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Abstract 

We consider a processor shared M/M/l queue that can accommo- 
date at most a finite number K of customers. We give an exact ex- 
pression for the sojourn time distribution in the finite capacity model, 
in terms of a Laplace transform. We then give the tail behavior, for 
the limit K — > oo, by locating the dominant singularity of the Laplace 
transform. 
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1 Introduction 

Processor sharing (PS) is one of the most interesting service disciplines in 
queueing theory. The M/M/l-PS queue assumes Poisson arrivals with rate 
A and exponential i.i.d. service times with mean 1///. The traffic intensity 
is p = X/fi. Here we consider a system which can accommodate at most K 
customers. If the system is filled to capacity, we assume that further arrivals 
are turned away and lost. 

Most past work on PS models deals with systems that have an infinite 
capacity of customers. In [2J, Coffman, Muntz, and Trotter analyzed the 
M/M/l-PS queue, and derived an expression for the Laplace transform 
of the sojourn time distribution, conditioned on both the number of other 
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customers seen, and the amount of service required, by an arriving customer. 
Using the results in [2], Morrison [6] obtains asymptotic expansions for the 
unconditional sojourn time distribution in the heavy traffic limit, where the 
Poisson arrival rate A is nearly equal to the service rate \x (thus p = X/fj, f 1). 

A service discipline seemingly unrelated to PS is that of random order 
service (ROS), where customers are chosen for service at random. In [7] 
Pollaczek derives an explicit integral representation for the generating func- 
tion of the waiting time distribution WroSj from which the following tail 
behavior is obtained 

Pr [W ROS > t] ~ e-^ 17 V 5A \ t -> oo. (1) 

Here a, @ and 7 are explicitly computed constants, with a = (1 — ^fp) 2 ■ 
Cohen [3] establishes the following relationship between the sojourn time 
Vps in the PS model and the waiting time Wros m the ROS model, 

pPr[V PS >t] =Pr[Wftos >t], (2) 

which extends also to the more general G/M/l case. In pQ relations of 
the form ([2]) are explored for other models, such as finite capacity queues, 
repairman problems, and networks. 

The present finite capacity model will have purely exponential behavior 
for t — > 00, and thus the subexponential (e~^* 1/3 ) and algebraic (t~ 5 / 6 ) 
factors that appear in (TTJ will be absent. Writing Pr[V > t] ~ C e~ St for 
t — > 00 for the finite capacity model, the relaxation rate 5 = 5{K) depends 
on the capacity, and we shall study its behavior for K — > 00 and for various 
values of p. In particular, if p < 1 and K — > 00 it will prove instructive to 
compare the tail behaviors when K = 00 (where ([1]), ([2]) apply) and when 
K is large but finite. 

Knessl [1] obtains asymptotic expansions for the first two conditional 
sojourn time moments for the finite capacity model. The assumption is that 
K S> 1 and separate analyses are carried out for the cases p < 1, p = 1 
and p > 1. In [5] he obtains expansions for the sojourn time distribution 
by performing asymptotic analyses for the three scales p — 1 = 0(iT -1 ); 
p - 1 = bK~ 1 / 2 , b > and p > 1. 

The main purpose of this note is to give an explicit exact expression for 
the Laplace transform of the sojourn time distribution in the finite capacity 
M/M/l-PS model. We apply a Laplace transform to the basic evolution 
equation and solve it using a discrete Green's function. The solution is 
summarized in Theorem 2.1. Then in Theorem 2.2 we locate the dominant 
singularity, which leads to the tail behavior, for K S> 1. 
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2 Summary of results 

We set the service rate p, = 1. Then the traffic intensity is p = A > 0. We 
define the conditional density of the sojourn time V by 

p n (t)dt = Pr [V G (t,t + dt)\N(0~) = n], < n < K - 1. 

Here N(0~) is the number of customers present in the system immediately 
before the tagged customer arrives. From [5], the quantity p n (t) satisfies the 
evolution equation 

p' n (t)=PPn+l(t)-(l + p)Pn(t) + — T Pn-l{t), < n < K - 1, (3) 

n + 1 

for t > 0, with the initial condition p n (0) = ^xj- If n = K — 1 the 
term pp n+ i(t) is absent, and we shall use the "artificial" boundary con- 
dition p K _ 1 {t) = p K (t). Taking the Laplace transform of (|3|) with p n {9) = 
Jo°° Pn(t)e~ et dt and multiplying by n + 1, we have 

(n + 1) pp n+1 (6)-(n + l) (l + p + e)p n (9) + np n - 1 {e) = -l (4) 

for n = 0, 1, K — 1 with the boundary condition 

P*-i(0) (5) 

Solving the recurrence equation with ([5]), we obtain the following 
result. 

Theorem 2.1 The Laplace- Stieltjes transform of the conditional sojourn 
time density has the following form: 

n K-i . „ K-1 

p n {9) = MG n Y P l Hi + MH n Y P l Gi - ~^ MH n Y ^ ( 6 ) 

/=0 l=n+l K 1=0 

for n = 0,1, K — 2 and 

F ft z=o 

Here 

M = M(9) = z^^y, (8) 
G n = G n {9) = z n {z + - z)- a (z^ - z) a - x dz, (9) 
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H n = H n {6) = — / z n (z + - z)- a {z - z-) a ~ x dz, (10) 

C is a closed contour in the complex z-plane that encircles the segment 
[z— , z+] of the real axis counterclockwise and 

AGk = Gk — Gk-i, 

AH K = H K - H K - X , (11) 
1 r 



z± = z±{9) = — l + p + e±y/(l + p + 9) 2 -4p , (12) 
2pL J 

a = a (9) = Z+ . (13) 

Z-Y- — Z— 

The singularities of p n (0) are poles, which are all real, and solutions of 
Hk{&) = Hk—i(0). A spectral expansion can be given in terms of these 
poles, but we believe that is more useful for certain asymptotic analyses, 
such as K — > oo. Below we give the least negative pole 9 S , which is also the 
tail exponent, i.e. limt^ooli -1 log Pr[V > t]}. 

Theorem 2.2 The dominant singularity 9 S in the 9 -plane has the following 
asymptotic expansions, for K — > oo ; 



1. P< 1, 



,i-Vp) ~it + k^~i5k^ + 0{k } ' (14) 



where ro ~ —2.3381 is the largest root of the Airy function Ai(z). 

2. p = l + r) K- 2 / 3 with r] = 0(1), 

lri 16r 3 + 8rj 2 r\ + (rf + I9r])n + r? 3 + 9 „,„_ 2 ^ 

9s -~K + K^ + ( } ' 

(15) 

where r\ = ri(rj) is the largest root of 

Ai'jn + ?? 2 /4) _ y 

Ai( ri + V 2 /A) 2' 1 ' 

3. p > 1, 

6s = -J<- (p-\)K 2 ~ (p - I)*** + (p-l)*K* + ° {K ~ 5) - (17) 
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Comparing e dst with (|14|) to ([T|) we see that they both contain the domi- 
nant factor e^^^v^) 2 * but whereas ([!]) has the sub exponential factor e - ^* 1 3 , 
(|14p leads to purely exponential correction terms that involve the maximal 
root of the Airy function. 



3 Brief derivations 

We use a discrete Green's function to derive © and ©. Consider the 
recurrence equation (jl]) and ©. The discrete Green's function Q(9;n,l) 
satisfies 

(n + l)pQ{9; n + 1, 1) - (n + 1)(1 + p + 9) G(9; n, I) 

+ n g(9;n- 1,1) = -S(n,l), (n, I = 0, 1, K - 1) (18) 

and 

G(9;K,l)=G(9;K-l,l), (I = 0,1, K), (19) 

where 5{n,l) = l{ n= /} is the Kronecker delta. To construct the Green's 
function requires that we have two linearly independent solutions to 

(n + l)pG(9;n + l,l)-{n + l)(l + p+9)G(9;n,l)+nG{9;n-l,l) =0, (20) 

which is the homogeneous version of (|18p . We seek solutions of (|20p in 
the form G n = J v z n g(z)dz, where the function g(z) and the path D of 
integration in the complex z-plane are to be determined. Using the above 
form in (|20p and integrating by parts yields 

z n g(z)[pz 2 -(l + p + 9)z+l)\ v 

- [ z n [{pz 2 -(l + p + 9)z + l)g'(z) + pzg(z)] dz = 0. (21) 

The first term represents contributions from the endpoints of the contour 
V. 

If (|21|) is to hold for all n = 0, 1, K — 1 the integrand must vanish, so 
that g{z) must satisfy the differential equation 

[pz 2 ~{l + p + 9)z + l] g\z) + pzg(z) = 0. 

We denote the roots of pz 2 — (1 + p + 9)z + 1 = by z + and with 
z + > Z- > for real 9. These are given by (|12p and if a is defined by 
(I13p . the solution for g(z) is g(z) = (z + — z)~ a (z_ — z) a_1 . By appropriate 
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choices of the contour T>, we obtain the two linearly independent solutions 
G n (cf. ([9])) and H n (cf. ([TO]) ), We note that G n decays as n — ► oo, and is 
asymptotically given by 

G„~^^ + "(z+-z_)- Q , rwoo. (22) 

However, G n becomes infinite as n — ► —1, and nG n -\ goes to a nonzero 
limit as n — > 0. Thus G n is not an acceptable solution at n = 0. -ff n is finite 
as n — ► —1, but grows as n — > oo, with 

F n ~^<+ 1 - a (z + -z_r- 1 , rwoo. (23) 
r(a) 

Thus, the discrete Green's function can be represented by 
A / AH n + BG n iil<n<K 

Here A, B and C depend upon 6, I and K. Using continuity of Q at n = I 
and the boundary condition (|19p . we obtain from (|24p 



HiAGk c b = ^Aff A - c 



HiAGk-GiAH k ' HiAGk-GiAH k 
Hence, (|24p can be rewritten as 

' H u AGk — G u AHk 



G(e-,n,i) 



CHr, if < n < / < K. 



To determine C, we let n = / in (|18|) and use the fact that both Gi and 
-ff; satisfy (|20p with n = I, which gives 

q GiAHk — HiAGk 

~~ (l + l) P AH K (G l Hi +1 -G l+1 H l y ( ' 

From (|20p we can infer a simple difference equation for the discrete Wron- 
skian Gi Hj+i — G7+1 whose solution we write as 

Gi H l+i - Gi + i Hi = (26) 
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where C\ = C\{9) depends upon 9 only. Letting / —* oo in (|26p and using 
the asymptotic results in (|22p and (|23p . we determine C\ as Ci = pM and 
then by (f25j) obtain 

C = — - -p l M, 

AH K 

where M = M(9) is defined by p. 

Then, we multiply (|18p by the solution pz(#) to ([5]) and sum over < 
I < K — 1. After some manipulation and using (fTUl) . this yields 

p„w = X)^ n »o. n = o,i,...,K-i, 

which is equivalent to © and ([7]). 

In the remainder of this section, we obtain the dominant singularity of 
p n (0) in the #-plane. From (0), the dominant singularity comes from solving 
AH K = H K - Hk^ = 0. 

We first consider p < 1. We evaluate H n in (|10p by branch cut integra- 
tion, which yields 

s.=^rrf-.-r'( !t -!)-c (27) 



sin(Q7r) f z+ £ — 1 / 6 — z_ 



Thus, we have 

Since the factor sin(a7r) appears in both the numerators and denominators in 
([6]) and (|7Dj zeros of sin(a7r) do not correspond to poles o£p n (9). To be at a 
pole we must have the integral in (|28p vanish. When p < 1 and the 
effects of finite capacity should become asymptotically negligible and the 
finite capacity model may be approximated by the corresponding infinite 
capacity model. Thus, from (pQ) the dominant singularity 9 S should be close 
to -a = -(1 - Jp) 2 . We scale 9 = -(1 - Jp) 2 + s/K and let 

i= z_ ± ±z. z_ ± -z_ w 

s 2 2 v ; 

in (|28p . After some calculation, the integral in (|28p is approximately given 
by 



^(w, s) + OiK- 1 ' 2 )} e v^/K*) dw 5 (30) 
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where 



ti \ -1/4 r , P 1/4 i ( l + w 
J(w,s)=p 1 w \/s + -^—^log 



g{w,s) 



2^T~ & Vl-W 
s(l-w 2 ) 



exp 



By solving -J^f(w,s) = 0, we find that there are two saddle points, at 
itf* = iy 7 ! + yfp/s. In order to be in the oscillatory range of (j30j) the two 
saddles must coalesce, which occurs at w = when s = —J~p. Then we 
introduce the scaling w = uK^ 1 ^ and s = y/p(—l + rK~ 1 / 3 ) and expand 
the functions f(w,s) and g(w,s) as 



f(w,s) 



i ru+- 



u 



9 q 

r u ru u 



K- 1 ' 2 -i { '-f + i|- + ^ ) i^ 5 / 6 +0 (E^ 7 / 6 ) , (31) 



SK s) = e" 1 ^ [l + (~ + u 2 )^ 1 ' 3 + 0(K- 2 / 3 ) 



(32) 



By using the definition of the Airy function 

1 f 00 

Ai( z ) = ± / e-W^dt, 
along with (|3ip and (|32p . the integral in (13QH can be evaluated as 



2vre" 1/2 K- 1 ' 6 \ AiM + K' 1 / 3 ± r Ai(r) + —r 2 Ai'{r) + 0{R- 2 / 3s 

L 15 15 



(33) 



To find the dominant singularity, (|33p must vanish. It follows that r is 
asymptotic to the maximal root of Ai. To get a refined approximation, we 
expand r as r = ro + aoK^ 1 / 3 + 0(K~ 2 / 3 ) and calculate cxq from ([33]) . Then 
using = -(1 - ^/p) 2 - ^/p/if + ^fprK^I 3 we obtain (fllj) . 

Next, we consider p = 1 + r\K~ 2 l 3 with 77 = 0(1). We now scale = 
- l/fif + R K~ 4 / 3 andw = « if" 1 / 6 , with R and u = O ( 1 ) . Then the leading 
order approximation to the integrand in (128ft is 



-_ftT 1 / 3 r ) /2-l/2^-5/6 



III- 



)+0{K- 1 / 3 )] e ^{-i[^-+{R+^)u]}. (34) 



Integrating (f34"j) over u now leads to 

2jre -K^r,/2-l/2 K -5/6 



\M{R + - + + 0(^-V3)] . (35) 



S 



Setting (|35p equal to zero we obtain the expansion R = r% + ai-fT -1 / 3 + 
0(-fT~ 2 / 3 ), which leads to # s ~ —1/K + riK~^ 3 , where r\ satisfies (fl~6j) . 
We note that to obtain the 0(K~ 5 ^ 3 ) term in (fT5|) . we needed the 0(iT -1 / 3 ) 
correction terms in f|34[) and (|35p . which we do not give here. 
Finally, we consider p > 1. We use (fTU|) and rewrite AiJ^- as 



A' 



2vri 



a-1 



(36) 



Scaling z = z+ + tt/i^, the contour C may be approximated by the contour 
£, which starts at — oo — iO, below the real axis, circles the origin in the 
counterclockwise direction and ends at — oo + i0, above the real axis. Thus, 
becomes 



AH 



K 



K c 



-1K-1 



\a-l 



h + h/K + h/K 2 + 0(K~ 



(37) 



where the are expressible in terms of the V function, with 



h = 



1 



u/z+ 



(z + — l)u a du 



z^{z + - 1)/I», 



and I 2 and h can be given as similar contour integrals. To find the dominant 
singularity, (|37ft must vanish. Viewing the Ij = Ij(9) as functions of 9, I\ 
will vanish when 24. = z+{9) = 1, which occurs at 9 = if p > 1. Thus we 
expand 9 as 6 = a 2 R- x + /3 2 R- 2 + 0{R- 3 ). Then §7$ yields 



fi(O)0 + \l" (0)9 2 + K- 1 [I 2 (0) + I' 2 (0)8] + K- 2 I 3 (0) = 0(K- 



so that a 2 



-/ 2 (0)//[(0) and 



02 



2 / i(0)«2+^(0)«2+/3(0)J// 1 (0), 



which leads to 9 S ~ — 1/K — l/[(p — 1)-K" 2 ]- We note that to obtain the 
higher order approximations in fjlT|> . we needed to compute two more terms 
in (1371). This concludes our derivation. 



References 

[1] S. C. Borst, O. J. Boxma, J. A. Morrison, and R. Nunez Queija, The 
equivalence between processor sharing and service in random order, 
Oper. Res. Letters 31 (2003), 254-262. 



9 



[2] E. G. Coffman, Jr., R. R. Muntz, and H. Trotter, Waiting time distri- 
butions for processor-sharing systems, J. ACM 17 (1970), 123-130. 

[3] J. W. Cohen, On processor sharing and random service (Letter to the 
editor), J. Appl. Prob. 21 (1984), 937. 

[4] C. Knessl, On finite capacity processor-shared queues, SIAM J. Appl. 
Math. 50 (1990), 264-287. 

[5] C. Knessl, On the sojourn time distribution in a finite capacity proces- 
sor shared queue, J. ACM 40 (1993), 1238-1301. 

[6] J. A. Morrison, Response-time distribution for a processor-sharing sys- 
tem, SIAM J. Appl. Math. 45 (1985), 152-167. 

[7] F. Pollaczek, La loi d'attente des appels telephoniques, C. R. Acad. Sci. 
Paris 222 (1946), 353-355. 



10 



